class: center, middle, inverse, title-slide .title[ # VEIN model and future integration with MUSICA ] .author[ ### Sergio Ibarra-Espinosa¹ ² ] .author[ ### 1 CIRES, University of Colorado-Boulder ] .author[ ### 2 GML, National Oceanic and Atmospheric Administration (NOAA) ] .date[ ### 2022-10-05 ] --- class: inverse center middle ## VEIN model and future integration with MUSICA > Sergio Ibarra-Espinosa¹'² <br> > 1 CIRES, University of Colorado-Boulder<br> > 2 GML, National Oceanic and Atmospheric Administration (NOAA)<br> > sergio.ibarra-espinosa@noaa.gov<br> > Presentation for MUSICA group, ACOM/NCAR <img src="https://atmoschem.github.io/vein/reference/figures/logo.png" alt="drawing" height="100"/> <img src="logo.png" alt="drawing" height="100"/> <img src="https://atmoschem.github.io/eixport/reference/figures/logo.gif" alt="drawing" height="100"/> <img src="https://www2.acom.ucar.edu/sites/default/files/styles/extra_large/public/images/MUSICAlogo_0.png" alt="drawing" height="100"/> <img src="https://ruc.noaa.gov/wrf/wrf-chem/wrfchem_logo.png" alt="drawing" height="100"/> --- # Intro .pull-left[ ## Why do we need MUSICA? Empirical and modeling studies have provided strong evidence of dynamical and chemical coupling across the range of spatial and temporal scales inherent in the Earth System. Current models, however, are inadequate in accounting for the two-way coupling of the range of relevant scales **from street-level to regional to global** and from the surface up to the mesosphere. ] .pull-right[ ## What is VEIN? [](http://cran.r-project.org/web/packages/vein) [](http://cran.r-project.org/package=vein) - R package to calculate vehicular emissions - Includes Fortran subroutines with // OpenMP - Detailed speciation applying Carter [(2015)](https://www.tandfonline.com/doi/full/10.1080/10962247.2015.1013646) - 13 direct publication [RG][https://www.researchgate.net/project/VEIN-An-R-package-for-vehicular-emissions-inventories] ] --- # vein  --- # Characteristics - Projects designed to produce comprehensiva inventories for people without knowledge in R - Imports sf package with bindings to GEOS, GDAL, PROJ and UDUNITS - Imports data.table: [super fast] https://h2oai.github.io/db-benchmark/ - main paper: https://gmd.copernicus.org/articles/11/2209/2018/ - 33K dl, 3732 article views, 28 crossref - YouTube Channel https://www.youtube.com/channel/UC2oYaS9mpnIDk8w55O8_bTg ### Process - Top down (POLYGON) and bottom-up (LINESTRING) - Activity adjusted by fuel consumption - Activity and emission factors by age of use - Includes methodolologies for China, US/EPA (MOVES), Brazil and Europe. - **Gridded emissions with any polygon shape** - Includes more than 10 chemical mechanisms --- # Atmoschem R packages .pull-left[ - **eixport** reads spatial emissions and generates inptus for air quality models - such as WRF CHEM - **respeciate** Port the US/EPA Speciate 5.0 to R (no need of MS Access) - **cptcity** (7120 colour palettes from cptcity web archive) ] .pull-right[ <img src="https://atmoschem.github.io/eixport/reference/figures/logo.gif" alt="drawing" width="100"/> [](http://cran.r-project.org/web/packages/eixport) [](http://cran.r-project.org/package=eixport) [](https://github.com/atmoschem/eixport) <img src="logo.png" alt="drawing" width="100"/> <img src="https://camo.githubusercontent.com/43369365a280712734504191f0cedc3a6fb34844686089b30c726179810704eb/68747470733a2f2f692e696d6775722e636f6d2f347639574d51762e676966" alt="drawing" width="100"/> [](http://cran.r-project.org/web/packages/cptcity) [](http://cran.r-project.org/package=cptcity) [](https://github.com/ibarraespinosa/cptcity) ] --- class: inverse center middle # Some applications --- class: center middle ## Formaldehyde: 0.059 Mt/y <img src="figuras/PM25.png" alt="drawing" height="550"/> --- class: center middle ## NP 50 nm: 1.58E+26 /y <img src="figuras/N_50nm_urban.png" alt="drawing" height="550"/> --- class: center middle ## NorthEast China <img src="figuras/fig7.png" alt="drawing" height="550"/> --- ## Applications: Sao Paulo .pull-left[ ```r library(sf) library(vein) x <- readRDS("rds/CO.rds") g <- st_sf( geometry = st_make_grid( x = x, cellsize = 3000, square = F)) co <- emis_grid(spobj = x, g = g) ``` ``` ## Sum of street emissions 1080803725.33 ## Sum of gridded emissions 1080803725.33 ``` ] .pull-right[ ```r plot(co[as.numeric(co$V9) > 0, "V9"], axes = T, lty = 0) ``` <img src="index_files/figure-html/unnamed-chunk-2-1.png" width="100%" style="display: block; margin: auto;" /> ] --- class: center ## WRF Chemi using eixport <video width="520" height="440" controls> <source src="wrfc.mp4" type="video/mp4" height="550"/> Your browser does not support the video tag. </video> --- class: inverse center middle # Integration with US/EPA MOVES --- ## Recent Research: Integration of MOVES and VEIN: - MOVES is the official vehicular emissions model for US. Runs on Windows, written in Java/SQL with MariaDB. VEIN is very versatile, ideal for traffic flow at streets. Currently has two approaches: - **1** Estimation using Windows with MOVES >3.0 installed. Emission factors are accessed using SQL in R. - **2** Estimation using any OS. Emission factors are exported from Windows as .csv.gz and read with `data.table::fread`. - Paper will be submitted to GMD (under development) --- ## Screenshots ```r vein::get_project(directory = "sacramento", case = "moves") ``` .pull-left[ <img src="figuras/main.png" width="83%" style="display: block; margin: auto;" /> ] .pull-right[ <img src="figuras/moves1.png" width="75%" style="display: block; margin: auto;" /> ] --- ## Sacramento County 2017 .left-column[ - Traffic flow for a 4-stage travel demand model output for Sacramento Area, extracted for Sacramento County.' - Traffic flow is for 2016 from CARB. - Traffic flow is total traffic volume 08:00-09:00. - Vehicular composition based on baltimore. - Fuel consumption for 2017. - Emission factors from Baltimore 2017. - Temporal factors from hourly VMT MOVES Baltimore. - Assumed BPR parameters. ] .right-column[
] --- class: center ## Speed parameters .pull-left[ <img src="figuras/SPEED.png" alt="drawing" height="500"/> ] .pull-right[ <img src="figuras/SPEEDBIN.png" alt="drawing" height="500"/> ] --- ## Emissions .pull-left[ <img src="figuras/emisac.png" width="93%" style="display: block; margin: auto;" /> ] .pull-right[ <img src="figuras/no2mov.png" width="93%" style="display: block; margin: auto;" /> ] --- ## Spatial Emissions ``` ## Sum of street emissions 10466599.96 ## Sum of gridded emissions 10466599.96 ```